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Abstract 


We intend to study a modified version of the planar Circular Restricted Three-Body Problem (CRTBP) by 
incorporating several perturbing parameters. We consider the bigger primary as an oblate spheroid and emitting 
radiation while the small primary has an elongated body. We also consider the perturbation from a disk-like 
structure encompassing this three-body system. First, we develop a mathematical model of this modified CRTBP. 
We have found there exist five equilibrium points in this modified CRTBP model, where three of them are collinear 
and the other two are non-collinear. Second, we apply our modified CRTBP model to the Sun-Haumea system by 
considering several values of each perturbing parameter. Through our numerical investigation, we have discovered 
that the incorporation of perturbing parameters has resulted in a shift in the equilibrium point positions of the Sun- 
Haumea system compared to their positions in the classical CRTBP. The stability of equilibrium points is 
investigated. We have shown that the collinear equilibrium points are unstable and the stability of non-collinear 
equilibrium points depends on the mass parameter и of the system. Unlike the classical case, non-collinear 
equilibrium points have both a maximum and minimum limit of u for achieving stability. We remark that the 
stability range of и in non-collinear equilibrium points depends on the perturbing parameters. In the context of the 
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Studying the Equilibrium Points of the Modified Circular Restricted Three- 


Sun-Haumea system, we have found that the non-collinear equilibrium points are stable. 
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1. Introduction 


Celestial mechanics plays an important role in understanding 
the dynamics of solar system bodies (see, e.g., Murray & 
Dermott 1999; Souchay & Dvorak 2010; Lei 2021; Pan & 
Hou 2022). One of the problems in celestial mechanics is the 
Circular Restricted Three-Body Problem (CRTBP). The study 
of CRTBP has the aim to investigate the movement of an 
infinitesimal object under the gravitational influence of two 
primaries that have a circular orbit around their center of mass. 
CRTBP has several applications, such as for deep space 
exploration and satellite navigation. The classical version of 
CRTBP assumes the primaries are point masses and it only 
considers the gravitational interaction between them. There are 
five equilibrium points in the case of planar. Three of them are 
collinear (Lı, L2, and L3) and the other two are поп-соШпеаг 
(L4 and Ls) (Murray & Dermott 1999). In order to make the 
CRTBP model more realistic, the classical version has been 
modified by considering several additional parameters. 

A stellar object, including the Sun, emits radiation. This 
radiation exerts pressure on objects in its path. There have 
been numerous studies that have considered radiation pressure 
force as another additional force in the restricted three-body 
problem (see, e.g., Haque & Ishwar 1995; Ishwar & Elipe 
2001; Kushvah et al. 2007; Kushvah 2008a; Das et al. 2009; 


planets and satellites: dynamical evolution and stability 


Yousuf & Kishor 2019; Patel et al. 2023). For instance, the first 
study on this topic was done by Radzievskii (1950). Chernikov 
(1970) extended the study by considering the relativistic 
Poynting-Robertson effect. Simmons et al. (1985) examined 
the effect of radiation pressure force in all ranges of value. 
More recently, Idrisi (2017) and Idrisi & Ullah (2018) 
considered the effect of planetary albedo on CRTBP as a 
consequence of solar radiation pressure force. 

Since the stars and planets are not perfectly spherical, 
another aspect that has been considered in the CRTBP is the 
oblateness of the primaries. Early studies about the impact of 
an oblate primary on the dynamics of restricted three-body 
problem were published by Danby (1965), Sharma & Subba 
Rao (1978, 1986). More recently, the effect of oblateness on 
the dynamics of CRTBP has been studied in detail by several 
authors (see, e.g., Markellos et al. 1996; Douskos & 
Markellos 2006; Safiya Beevi & Sharma 2012; Abouelmagd 
et al. 2013; Zotos 2015; Yousuf et al. 2022). Moreover, some 
authors have considered the effect of both oblateness and 
radiation force in their calculation. For instance, Singh & 
Ishwar (1999) studied the linear stability of triangular 
equilibrium points when both primaries are oblate and emitting 
radiation. This study was extended by Singh (2009) for the 
nonlinear stability of L4. AbdulRaheem & Singh (2006) 
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Table 1 
The Abscissa Position of Collinear Equilibrium Points (L4, 15, and L4) in Sun-Haumea System with и = 2 x 10^? and T = 0.11 

1— qd A 1 M, Lı 1» La 

1 0 0 0 —1.000873832771965 —0.999126671989864 1.000000000833333 
1.6 x 1075 2.6 x 107!! 3.5 x 1077 3x107 —1.00087355691303 —0.999126395886954 0.999999369260791 
1.6 x 1074 2.6 x 107!! 3.5 x 1077 3x107 —1.00085632697634 —0.999108413962744 0.999946566436973 
1.6 x 107? 2.6 x 107!! 3.5 x 1077 3х 107' —1.00087373433354 —0.999126573666542 0.999999902060868 
1.6 x 107 2.6 x 107!! 3.5 x 1077 3 x 107' —1.00087355691303 —0.999126395886954 0.999999369260791 
1.6 x 107 2.6 x 107? 3.5 x 1077 3x107 —1.00087355691116 —0.999126395888830 0.999999369260793 
1.6 x 107 2.6 x 10713 3.5 x 1077 3x107 —1.00087355691305 —0.999126395886936 0.999999369260791 
1.6 x 1079 2.6 x 1071! 3.5 x 1077 3 x 1077 —1.00087355691303 —0.999 126395886954 0.999999369260791 
1.6 x 107 2.6 x 107!! 3.5 x 10? 3х 107' —1.00087402445000 —0.999125928668234 0.999999368852499 
1.6 x 1075 2.6 x 107!! 3.5 x 10? 3x107 —1.00087355686628 —0.999126395933676 0.999999369260832 
1.6 x 107 2.6 x 107!! 3.5 x 1077 3х 107' —1.00087355691303 —0.999126395886954 0.999999369260791 
1.6 x 107 2.6 x 107!! 3.5 x 1077 3 x 10? —1.00086393713086 —0.999116570293592 0.999989644085257 
1.6 x 1075 2.6 x 107!! 3.5 x 1077 3x 10° —1.00087365419779 —0.999126493044322 0.999999466517286 


investigated the dynamics of CRTBP when both of the 
primaries are oblate and emit radiation, together with 
perturbations in the Coriolis and centrifugal force. Other 
authors, such as Nurul Huda et al. (2015), Dermawan et al. 
(2015) and Mia et al. (2023), considered the effect of 
oblateness and radiation force in the Elliptic Restricted 
Three-Body Problem. 

Our solar system contains several types of celestial bodies. 
Among them are elongated objects like a few asteroids, comets, 
and dwarf planets. These celestial bodies can be approximately 
described as finite straight segments. Previous studies of 
CRTBP have been enriched by assuming one or both primaries 
have an elongated body. At first, Riaguas et al. (1999) and 
Riaguas et al. (2001) analyzed the dynamics of a two-body 
problem by considering one of the primaries as a finite straight 
segment. These works were extended by, e.g., Jain & Sinha 
(2014), Kaur et al. (2020), and Kumar et al. (2019), into the 
restricted three body-problem assuming both or one of the 
primaries have elongated shapes. In more recent studies, Verma 
et al. (2023a) examined the perturbed restricted three-body 
problem, where the smaller primary has an elongated shape and 
the larger primary is oblate and emits radiation. Verma et al. 
(2023b) considered the effect of a finite straight segment and 
oblateness to study the dynamics of the restricted 2 + 2 body 
problem. 

Meanwhile, the effect of a disk-like structure as a perturbing 
force near a three-body system has been well studied by several 
authors (see, e.g., Jiang & Yeh 2004; Kushvah 2008b; 
Kushvah et al. 2012; Kishor & Kushvah 2013; Mahato et al. 
2022a). Jiang & Yeh (2004) considered CRTBP by analyzing 
the influence of a disk-like structure near the three-body 
system. Yousuf & Kishor (2019) analyzed the effect of a disk- 
like structure, oblateness, and albedo on the CRTBP. Mahato 
et al. (2022a) extended the study of classical CRTBP by 


considering a disk-like structure and an elongated body. 
Mahato et al. (2022b) investigated the stability of equilibrium 
points within a framework of the perturbed restricted 2 + 2 
body problem, taking into account the influence of a disk-like 
structure. 

This study aims to obtain the collinear and non-collinear 
equilibrium points and investigate their stability under a 
framework of modified CRTBP incorporating the effect of 
radiation pressure, oblateness, finite straight segment, and disk- 
like structure. We intend to extend the work of Yousuf & 
Kishor (2019) by assuming the small primary to be a finite 
straight segment rather than being oblate. It is also an extension 
of Mahato et al. (2022a) since we consider the effect of 
oblateness and radiation from the bigger primary. 

Here we apply our modified CRTBP model to the Sun— 
Haumea system by assuming the Sun is a bigger primary with 
an oblate shape and emitting radiation and Haumea is a smaller 
primary which has an elongated body. We also consider the 
Kuiper Belt as a disk-like structure surrounding the Sun- 
Haumea system. Haumea was chosen as our case study because 
of its unique characteristics, which have captured the attention 
of scientists since its discovery in 2003. The surface of Haumea 
is dominantly covered by water ice (Barkume et al. 2006; 
Pinilla-Alonso et al. 2009; Noviello et al. 2022). There is also 
evidence that organic material exists оп Haumea's surface 
(Lacerda et al. 2008; Gourgeot et al. 2016). Recently, it was 
discovered that Haumea has a ring and two satellites named 
Namaka and Hi'iaka (Ortiz et al. 2017). Moreover, previous 
studies have proposed Haumea as a destination for space 
missions in the coming decades (see, e.g., Grundy et al. 2009; 
Sanchez et al. 2014). 

Besides the Sun-Haumea system, this modified CRTBP 
model can be applied to other cases. For instance, many 
planetary systems outside of our solar system have been 
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Table 2 
7? and T= 0.11 
1-4 A І M, 145 
x y 

1 0 0 0 —0.499999998000000 +0.866025403784439 
1.6 x 1079 2.6 x 107 !! 3.5 x 1077 3 x 1077 —0.499999464678626 +0.866024982450545 
1.6 x 10-4 2.6 x 107"! 3.5 x 1077 3x 10-7 —0.499946656129219 +0.865994492868782 
1.6 x 10° 2.6 x 107 !! 3.5 x 1077 3х 10-7 —0.499999997479636 +0.866025290063446 
1.6 x 107 2.6 x 1071! 3.5 x 1077 3x 107 —0.499999464678626 +0.866024982450545 
1.6 x 1079 2.6 x 107? 3.5 x 1077 3 x 1077 —0.499999465965624 +0.86602498 1707493 
1.6 x 107 2.6 x 10713 3.5 x 1077 3х 107' —0.499999464665756 +0.866024982457975 
1.6 x 107 2.6 x 107"! 3.5 x 1077 3x 107 —0.499999464678626 +0.866024982450545 
1.6 x 107 2.6 x 107"! 3.5 x 10? 3 x 1077 —0.499999464372405 +0.866024982155884 
1.6 x 1079 2.6 x 1071! 3.5 x 10? 3 x 1077 —0.499999464678656 +0.866024982450574 
1.6 x 107 2.6 x 107"! 3.5 x 1077 3 x1077 —0.499999464678626 +0.866024982450545 
1.6 x 107 2.6 x 107!! 3.5 x 1077 3x 107 —0.499999464663069 +0.866013755215885 
1.6 x 107 2.6 x 107"! 3.5 x 1077 3x 10° —0.499999464678781 +0.866025094722156 


discovered, and some systems have been found to have dust 
particle disks or asteroid belts, which are believed to be similar 
to the Kuiper Belt or main belt in our solar system (see, е.2., 
Greaves et al. 1998; Matra et al. 2019). Meanwhile, previous 
studies explained the presence of extrasolar asteroids or dwarf 
planets near the host star (see, e.g., Jura 2003; Dufour et al. 
2010). Moreover, some space explorations have been devoted 
to exploring small solar system bodies near the main belt or 
Kuiper Belt region. It is known that several solar system bodies 
have an irregular shape. Therefore, it is reasonable to study the 
combined effects of perturbations from a disk, an elongated 
body, and an oblate radiating body on the motion of an 
infinitesimal mass in the CRTBP. 

The structure of this paper is as follows. In the next section, 
we present a mathematical formulation of the dynamical model. 
The position and stability of equilibrium points are elucidated 
in Section 3. Section 4 describes the implementation of the 
dynamical model in the Sun—Haumea system. Finally, the 
conclusion is provided in Section 5. Here, MATLAB’s 
Symbolic Toolbox is used to conduct certain algebraic 
calculations and find numerical solutions. 


2. Mathematical Formulation of the Dynamical 
System 


In this work, we consider a system where an infinitesimal 
mass moves under the influence of a bigger primary with mass 
m, and a small primary with mass т». The primaries of this 
system have circular orbits around their center of mass. We 
treat the bigger primary as a source of radiation with an oblate 
spheroid shape, while the small primary has an elongated 
shape. The unit of time is normalized to make the Gaussian 
constant of gravitation equal to one. The mass parameter is 


represented by u=m/(m +m) where m,—1- и and 
то = p. In the case of a restricted three-body problem, it is 
more convenient to introduce the system in the rotational 
coordinate Oxy. The primaries are located on the x-axis with the 
distance between primaries chosen as the unit of length. The 
coordinates of the bigger primary, small primary, and the third 
body аге (и, 0), (и— 1, 0), and (x, y), respectively. The 
oblateness factor of the bigger primary can be represented by 
A= (АЕ? — АР?) / 5R? where A < 1, AE and AP represent the 
equatorial and polar radii, respectively, and R is the effective 
radius when assuming the primary to be a spherical object. 
Meanwhile, the radiation force F, acts opposite to the 
gravitational force and diminishes with respect to distance. 
The total force acting on the bigger primary can be written as 
F,—F,=4qF,, hence q—1 —(F,/F,). Here q is called the 
mass reduction factor where 0 < 1 — q «& 1. The small primary 
is assumed to be a finite straight segment with length 2/. The 
effect of a disk-like structure surrounding the system is also 
considered in this study. Following Miyamoto & Nagai (1975), 
the planar version of unitless potential disk-like structure is 
given by V(x, y) — M,/ Vr? +, T^. where M, is the total 
mass of disk-like structure, 7? = x? + у is the radial distance of 
the infinitesimal mass, and T = a + b is the total of flatness and 
core parameters. Let the distance of primaries to the center of 
mass be s, and s2. Considering the previous works such as 
Kushvah (2008b), Yousuf & Kishor (2019), and Mahato et al. 
(2022a), the motion of the primaries is given by 


— Gmm 1 3A GMpmr- 
| 2. р 282) (02 + TP 
Gmm 3A GM,m»r. 
2__ тту, , 22 | | Ит 
msn = С" (1 | meu 0 
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where R=s,+52 is the distance between primaries. 
p = 1— и + 2 means a dimensionless quantity of the 
reference radius of the disk-like structure (Singh & Taura 2014). 
Assuming R= 1, G= 1, and m, + m»; = 1, the mean motion n 
of the system can be calculated by adding both equations in 
Equation (1), approximating the expression 1/(1 — Ê) in series 
as 1 + P, and neglecting the term AP. Hence we have 
3 2Мьғ, 


па = 1 + 1 + ДА + 


2 (+P M 


Equations of motion of the third object in CRTBP are stated as 
follows 


00 
"End 
X ny Әх 
pyre (3) 
ду 


where U is a pseudo-potential function 


2 1— 1— pA 
п лед + у) + 400 + Copa 


ri ћ 
H тр + тэл + 21 Mp 
| lo | . 4 
21 ЈЕ + mm Jr? + T? e 


Here r} = x — u4-1—0? + y andr, = (х = gu + 1 + 52 + y? 
are the distance of a third body to the small primary and 
r? = (х — и)? + y? is the distance between the third body and 
the bigger primary. It should be noted that the equation of 
motion differs from the equation of motion in Yousuf & Kishor 
(2019) since, in our case, we regard the small body as a finite 
straight segment. 


3. Equilibrium Points 
3.1. Position of Equilibrium Points 


The conditions of equilibrium points are = y —X—y-— 
0. Hence we can deduce that Q, = Q, = 0, i.e., 


2, _ 44 — ша — 3(1 н) (х — WAG 


n 
n 2? 
2u х-ш+1- 1 х-и+1+1 
2 2 + 
(ra + r2) — 41 T2 T2 
Myx _ 
(r? ES T?)3/2 P 
(5) 
E qQ — ију _ 3 — p)yAg 
y 3 5 
r 2n 
2u y y Mpy = 
(па + r2 —4P\r r2 (r? + Т2)3/2 
(6) 
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In the following, we solve Equations (5) and (6) to find the 
position of equilibrium points. 

The collinear points are located in a line with the primaries, 
thus we have y — 0. Equation (5) becomes 


Q,(x, 0) = nix — ы ы; 
Ix — up 
— 30 — ш(х — ИЈА4 
2|x — up 


поа лиши пад ee E 
Ix-pti-Hün |х-и+1+/[ 


2 
=a 4 |> 0. 


х? 2x? 


(7) 


In order to find the solution, we divide the region into three 
parts, i.e., (C oo, w—1—J, (и— 1 — L и), and (p, оо). Here 
Li, 15, and L} are the solution located in (— co, u — 1 — D, 
(и — 1 — L и), and (u, оо), respectively. Hence we have 


me Meli AE (x = uy (x — 2и + 2)? — 4D) 
x 2x? 


+4@ — шб — uy (Ох – и + 2? — 4D) 
+З да — (Cx — Зи + 2° — 4) + их и), 


if -oo«x«pnu-—1-l 


[s Hfi з) p)*((2x — 2и + 2)? — 42) 
2x“ 


Ол. 0) = +901 — шб — uy (Qx – 2u + 2)? — 41) 
Зала = p)((2x — 2u + 2)? — 42) — 4u(x — uy, 


Кџи—1-+1 << х< и 
2 
(s al T )e p)*((2x — 2и + 2)? — AP) 
{2 


x 


q( — we — u)KQx — 2u + 2): — 40) 


3 
-344ü = p)((2x — 2u + 2)? — 4P) — 4u(x — uy. 
if y < x < oo. 


(8) 


These three equations have been solved numerically to find 
each collinear equilibrium point. Only the real solution is 
considered for the position of equilibrium points. 

Meanwhile, there are two non-collinear equilibrium points, 
i.e., L4 and Ls. The additional condition of these equilibrium 
points is y+ 0. Equations (5) and (6) can be rewritten in the 
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form 


2  40— и) 30 — u)Aq ги 
ud c 3 5 2_ APR 
ћ 2n (r21 + r22) 4l 


Mpx jus — н) | Зна = Ад 


(r? + T2572? " 2 
2u и+1—1 и+1-+1 
J 5 | = 0. 
(пр + 722)" — 41 T21 Г22 
(9) 
2  41— и) _ 3 — u)Aq 2u 
He =a t 2 2 
ri 2п (r21 + г2) — 4l 
x |r е 
721 722 
Му " 
(r? + Т2)3/2 
(10) 
Hence from Equation (10) we have 
„2 id= _ За – Ад ___ ги ___ 
n 2n (ro, + ry — 41 
[Se E 
Гл 722 
Му E 
(r? + Т2)3/2 
(11) 
Substituting Equation (11) into Equation (9) gives 
441 – и) , 30 — u)Aq 2 
И dw (га + ry? – 4P 
«( pal = nti tt) (12) 
721 722 


In the classical case, the positions of these equilibrium points 
are located at r; = 1 and rz = 1. Since some perturbations exist, 
we assume that r, and љ are perturbed by e; and е 
respectively. Hence, in our case, we have (Mahato et al. 2022a) 

r = 1+ €1; 


тоу = l+ & — 1/2; rn = l + & + 1/2. 


(13) 


The calculations of є and є are done by substituting 
Equation (13) into Equation (12) and Equation (11) and 
solving these equations. By approximating with series and 
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neglecting higher order terms of є, є, Ё, and A, we have: 


qom cR 24a e SE ҖЕ +2Ан 
4 q (p — 1) 
Syu- P (2 414) 
а (132 – 12) (u - 1) ' 
e= 16п - 40P y + 280 — 16 
52 P — 48 


where y = 1 + 12 + ЗА/2 + M,Qr, — 1)/(r? + T?)3/2. The 
position of non-collinear equilibrium points (x,, Yo) is given by 


3 
у= + g Be (15) 


Putting the value of e; > into Equation (15), we get: 


(14) 


= 2 2 _ 
wap 4 167- 4012 у + 2812 – 16 


2 52 12 — 48 
4+ 4 y LL 4u 
= – = –— 2А ад +744 4% +2A ра 
4q (и = 1) 


2(2и | Пор 
5yu-—L e p ) 
qı (1312 – 12) (а — 1) | 
Е 16у — 4012 y+ 2812 — 16 


ima | 52 12 — 48 
на FE уде doti SERA 
Е 3 14 3 3 1H 4 
4q,(u— 1) 


2(2p | Птн 1/2 
Ed 2 Па ) 
а 1312 – 12) (2 – 1) | ` 
(16) 


If the perturbation parameters are not considered, Equation (16) 
is similar to the classical version where x, = р — Н апа 


13 
у 


3.2. Linear Stability 
Let us assume a small displacement in an equilibrium point 
by defining 


и=х— Xo; y—y-y, (17) 


where “о” corresponds to the equilibrium points. The equation 
of motion from this small displacement is expressed as follows: 


ü — 2пў=иО + Убу, 
B+ 2пи = и ду, + УОу (18) 


уу, 
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where 


= = 2 — 
2:304 За ше ш ql : p) 
r т 
15А4(1 — )x = ш)? _ За — WA , 3Mpx? 
2r] 2 (?+т?у? 


M, 
(r? zi Т2)3/2 


| 2n 
(тә + r5)? — 4 


1 1 
x ——— O E 
[.——; у 


21 722 
Qu i, 4 
(га + тә)? — al от» 
(= | |, 
rj rip | 
(19) 
оо = 12 + 30 = шу? _ 20 — и) , 15А4(1— py? 
yy М ү 277 
3(1 — Ад 3M,y* M, 


P 2u | 1 1 ) 
(тә + r2) — APX ritro 1 та + r2 + 21 


Qu 1,1 у? у? 
(rai r3) — AP| гл Р 5 2%) ( 


оо _ 340 — WO — Wy , 15Aq(1 — WE — Wy 
xy 5 7 
ri 2n 
3Mpxy 
(г2 + “ 


1 
(г21 + — alas: + 722-21 ray + roc z) 


«(LL | l x pti 
721 722 721 722 


2 
+ ue - 
(ra + г2) — 41 


ДЕ +11 x zd 
rj 2» 


(21) 


Nurul Huda et al. 


Here 2° means the pseudo-potential is evaluated at equilibrium 
points. Hence it is constant. Equation (18) has general solutions 


4 
v= у) е^ (22) 


where o; and f; are constants while A; is the root of the 
characteristic equation. Substituting Equation (22) into 
Equation (18) produces 


Y= 0%, —2лА— 05, 
| Ој e 
2пх-0% X- РЕ Lo 


The first term on the left-hand side has to be a singular matrix. 
Hence the determinant of this matrix has to be zero 


X (4n? — 0%, – 05)% + 08,08 — (02): =0. (24) 
This equation is called the characteristic equation. It is a 
quadratic equation in А2. The solution of this quadratic 


equation is A;— Nam Jb? — 4c)/2, where b= 
4n? — (а, — OS, and c = Q2,05, — (0%)7. If all obtained 
A; are purely imaginary, then it; gives the motion of stable 
periodic behavior near the vicinity of equilibrium points. 
However, if there is at least one A; which has the form of real or 
complex, then the third body is unstable since u and v will 
exponentially increase with respect to time. We can investigate 
the stability behavior of the system by looking at the sign of b 
and c. The system is stable if b  O, b? — 4c > 0, and b > 


XP? — 4с since this case produces all pure imaginary А. 


4. The Case of Sun-Haumea System 


In this work, we model the Sun-Haumea system through the 
framework of the restricted three-body problem with the Sun as 
the bigger primary and Haumea as the small primary. Here we 
also consider the Kuiper Belt in this system. We assume 
Haumea has a circular orbit and orbits in the same plane as the 
Kuiper Belt. The mass of the small primary is a combination of 
Haumea's mass and the mass of Haumea's satellites: Namaka 
and Hi'iaka. The Sun has a mass around 1.989 x 10°° kg. 
Haumea has a length of ~2300 km for its largest axis and a 
mass of 4 x 10?! kg (Ragozzine & Brown 2009). Meanwhile, 
Namaka and Hi'iaka have masses of 1.79 x 10!? kg and 
17.9 x 10!? kg, respectively (Ortiz et al. 2017). Hence we have 
џи —2 x 107° and 1= 3.5 х 1077. Following Yousuf & Kishor 
(2019), here we assume that the Sun has A — 2.6 x 1071! while 
the Kuiper Belt has T= 0.11 and M, = 3 x 10 7. According to 
Sharma (1987), the photogravitational parameter q can be 
expressed in the CGS unit system as q = 1 — (5.6 x 107° /ар) 
where а апа р are the radius and density of a moving body, 


1—4 А 

1 0 

16 х 109. 26x10! 
1.6 x 1078 2.6х 107"! 
16х10 26x10! 
1.6 x 1076 26x10! 
1.6 x 1076 2.6х 107° 
1.6 x 1056 2.6 x 10773 
1.6 x 1076 26x10! 
1.6 x 1076 26x10! 
1.6 x 1076 26x10! 
1.6 x 1076 2.6 x 107"! 
1.6 x 1076 26x10! 
1.6 x 1076 2.6 x 107"! 


0 


3.5 x 
3.5 x 
3.5 x 


3.5 x 
3.5 x 
3.5 x 


3.5 x 
3.5 x 
3.5 x 


3.5 x 
3.5 x 
3.5 x 


1077 
1077 
1077 


1077 
1077 
1077 


1077 
107 
107° 


1077 
1077 
1077 


M, 


1077 
107 
1077 


1077 
1075 


Table 3 
Characteristic Roots of Collinear Equilibrium Points in Sun-Haumea System with и = 2 x 107° 
Li 1 Із 
А2 Азд А2 Азд А2 Азд 
2.50618628025287 2.070315207902671 2.51039039369088 2.072875380165811 0.000072456881366 1.00000000175000i 
2.50732613894453 . 2.07100940502565i 2.50924970700464 . 2.07218079713681i 0.000074119040096 1.00000030173167i 
2.58032763652133 2.115614498484641 2.43683616339914 2.028236060452701  0.000074121110154 1.000000301733141 
2.50659334481697 . 2.07056321377481i 2.50998449730882 . 2.07262831852643i 0.000074119018377 1.00000030173165i 
2.50732613894453 . 2.07100940502565i 2.50924970700464 . 2.07218079713681i 0.000074119040096 1.00000030173167i 
2.50732614822158 2.071009410679031 2.50924971630782 . 2.07218080279818; . 0.000074119037849 . 1.000000299801 167 
2.50732613885112 . 2.07100940496873i 2.50924970691157 . 2.07218079708017i 0.000074119039347 1.000000301750971 
2.50732613894453 . 2.07100940502565i 2.50924970700464 . 2.07218079713681i 0.000074119040096 1.00000030173167i 
2.50925730691060 2.072185425990871 2.51118123266234 2.073357254001331 0.000074119034105 1.00000030234410i 
2.50732594585502 2.071009287451131 2.50924951387951 . 2.07218067951878i 0.000074119040096 1.00000030173161; 
2.50732613894453 . 2.07100940502565i 2.50924970700464 4 2.07218079713681i 0.000074119040096 1.00000030173167i 
2.54765547950673 2.095629552586511 2.46924636836900 2.047881423101191  0.000172088738229 1.000030001421401 
2.50692401885402 . 2.07076439389322i 2.50965096408882 . 2.07242501730710i 0.000072473711215 1.000000004730571 


107? 


Note. Here i Means ¥—1. We used T = 0.11. 
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(c) L3 


Figure 1. Plot of и vs. characteristic roots (Ај 234) for L1, L2, and L3, with / = 3.5 x 1077, M; = 3 x 1077, A =2.6 x 10! , and 1 — а=16х 1079, The real and 
imaginary parts of characteristic roots are marked by solid and dotted lines, respectively. Here we used T — 0.11. 


respectively. Assuming a spacecraft has а = 700cm and 
p—0.05grcem ^, 1 — q = 1.6 x 1079. 

We calculated the position of the collinear equilibrium points 
of the Sun-Haumea system. By substituting the property of the 
system into Equation (8) and solving it numerically, we found 
Іл, Lọ, and L4. Table 1 shows the position of collinear 
equilibrium points. Here we vary the value of each perturbation 
parameter to examine the impact on the equilibrium point 
position. In the case of Lj, the position gets closer to the 
primaries if A and 1 — q increase. Decreasing A and increasing 
1 — q make Гл be closer to the bigger primary. The position of 
L4 is nearer with respect to primaries if the bigger primary 
emits stronger radiation pressure. According to Table 1, the 
position of collinear equilibrium points depends on the value of 
M, and I. Increasing M, and decreasing / make the location of 
Lı become nearer to the smaller primary. The increment of M, 
and / affects the position of L> to become closer to the bigger 


primary. L5 gets closer to the primaries if we increase the value 
of M,. 

The positions of non-collinear equilibrium points are 
calculated from Equation (16). Table 2 shows the positions 
of non-collinear equilibrium points with respect to the chosen 
value of several parameters. When there are no perturbing 
factors, the triangular points have the same coordinates as in 
the classical case. The inclusion of perturbation parameters has 
resulted in a shift in the location of non-collinear equilibrium 
points. The increment of A makes the position of these 
equilibrium points closer to the small primary. In contrast, if we 
reduce q or increase M,, the positions of equilibrium points are 
shifted toward the bigger primary. The positions are closer to 
the bigger primary in line with the increase of /. 

We now analyze the linear stability of each equilibrium point 
in the Sun-Haumea system. Collinear equilibrium points lie in 
the abscissa. Hence we have (2%, = 0. In order to study the 
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Figure 2. Plot of џ vs. characteristic roots (А) 23,4) in L4 with different parameter configurations. The real and imaginary parts of characteristic roots are marked by 
solid and dashed lines respectively. Here we use T = 0.11. The details of the parameters that are used in each subfigure are as follows. (a) A = 0, 1 — q = 0, 1= 0, 
M, = 0. (b) A=2.6 x 1071, 1 — q = 1.6 x 10°, 1= 3.5 x 1077, M, 23x1077. (c) A =2.6 x 10$ 1 — q = 1.6 x 10.5, L= 3.5 x 107, M, —3 x 107”. (d) 
А=2.6х10°!!, 1-4=1.6х 107, 1=35х 107, М=3 х 1077. (e) A=26x10'', 1-g=16x10°%, 1235x107, М=3х 107. (f) 


A=2.6 x 107"!, 1 — q = 1.6 x 107%, 123.5 x 10 5, M, = 3 x 1077. 
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stability, we divide the abscissa into three regions, i.e., Lı 
(œ, и—1—%, L (и—1—1, и), and Lz (и, оо), and 
calculate the sign of b and b? — 4c numerically for each region. 
First, we estimate the stability by considering the perturbation 
parameters in the Sun-Haumea system. As shown in Figure 1, 
there exist pure real and pure imaginary characteristic roots for 
и between 0 and 0.5. Hence, all collinear equilibrium points of 
the Sun-Haumea system are unstable. Furthermore, we 
conducted the calculation by varying the value of perturbation 
parameters. Table 3 displays the result of the calculation. All 
regions have b <0 and 22 — 4c > 0 which mean it produces 
two real pairs and two pure imaginary pairs. This shows that 
even if we change the value of perturbation parameters, the 
collinear equilibrium points remain unstable. 

Next, we investigate the stability of non-collinear equili- 
brium points in the Sun-Haumea system. We discuss only L4 
since the dynamics of Ls is nearly similar. In the classical case, 
non-collinear equilibrium points are stable under the condition 
27u(l— џи) <1. Hence we can deduce p< д, where the 
critical mass pu. = 0.038520896504551. This critical mass can 
be calculated by finding the solution of b? — 4c — 0. In this 
modified version of CRTBP, we numerically calculate the roots 
by solving Equation (24). By considering the perturbing 
parameters, it shows that the stability of non-collinear 
equilibrium points has a maximum limit (4) and minimum 
limit (4o) of mass parameters, an aspect which is different from 
the classical case. For the Sun-Haumea system, we found 
He = 0.0385208896007 and py, = 1.386 x 10-7, Since the 
Sun-Haumea system has и =2 x 107°, we conclude that the 
Sun-Haumea system has stable non-collinear equilibrium 
points. Figure 2 displays a comparison of stability for several 
cases by changing the perturbing parameters of the Sun- 
Haumea system. It shows that the range of stability depends on 
the parameters A, q, 1, and М,. The characteristic roots have the 
form of pure imaginary if j4 «i « He The considered 
perturbation parameters alter the range of stability in д. The 
increment of A or reduction of q reduces the size of the stability 
area. The stability region is shifted toward bigger u if M, and / 
increase. 


5. Conclusion 


We have investigated the dynamics of an infinitesimal mass 
under the gravitational influence of two primaries. Our study 
assumes that the smaller primary is an elongated body, while 
the larger primary is oblate and also emits radiation. In 
addition, we have taken into account the presence of a disk that 
surrounds the three-body system. We have found that there are 
five equilibrium points in this modified CRTBP such that three 
of them are collinear and the other two are non-collinear. Our 
numerical exploration of the Sun-Haumea system has revealed 
that the inclusion of perturbing parameters has caused a 
displacement in the position of the Sun-Haumea system's 
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equilibrium points with respect to their positions in the classical 
CRTBP. We noticed that the magnitude of the perturbing 
parameters (q, A, 1, and M,) can affect the positions of the five 
equilibrium points. This analysis shows that the non-collinear 
equilibrium points of the Sun-Haumea system are stable, while 
all collinear equilibrium points are unstable. Moreover, we 
have ascertained that the collinear equilibrium points remain 
unstable for several possible ranges of perturbing parameters. 
In contrast, the non-collinear equilibrium points are condition- 
ally stable with respect to u. When taking into account the 
perturbing parameters, we have found that there are upper and 
lower limits of и for achieving stability of non-collinear 
equilibrium points. The stability region for и depends on the 
perturbing parameters. 
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